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Abstract. 

The functional form of the nuclear symmetry energy has only been determined 
in a very narrow range of densities. Uncertainties concern both the low as well as the 
high density behaviour of this function. In this work different shapes of the symmetry 
^> ■ energy, consistent with the experimental data, were introduced and their consequences 

I for the crustal properties of neutron stars are presented. The resulting models are in 

agreement with astrophysical observations. 

O; 

1. introduction 

The energy per particle used to described neutron star interios can be expressed in terms 
of baryon number density n = np + nn and isospin asymmetry a - of the system: 

00 

!>■ ■ E{n, a) = V{n) + E,{n) + Oia") (1) 

o_ 

If instead of a the proton fraction x is introduced then then a = {l-2x), which proves to 
. be useful. In the model used here the only constituents of stellar matter are nucleons and 

^ ■ leptons: electrons and muons. Around and above the nuclear density n^ - 0.16 fm~^ 

nucleons and leptons form a quantum liquid, which stands for liquid core of the neutron 
star. Slightly below no matter cannot exist as a homogeneous fluid - the one-phase sys- 
tem is unstable and the coexistence of two phases is required. At these densities matter 
^ I clusterizes into positive nuclei immersed in a quasi-free gas of neutrons and electrons, 

most likely forming a Coulomb lattice with solid state properties and corresponds to 
the crust covering the liquid core of a star. 

For typical NS masses, between 1-2 M©, most of the stellar matter is occupied by 
the core, so the global parameters like the mass, radius, moment of inertia are com- 
pletely determined by the functional form of the Eq. ([T|l. Whereas the isoscalar part 
V{n) corresponds mainly for the stilfness of Equation of State (EOS) which is relevant 
for the maximum mass of NS, the isovector part Ef^jn) is respo nsible for the chemical 
composition of the matter (see (' Kubis & AIvarez-Castilloll2012h for details). Both func- 
tions V{n) and E An) have been irn plemented by means of Bezier functions composed 
of control points ( Wikipedial |20 1 2h 
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where ("^ is the binomial coefficient and the n + \ control points P,, (/ - 0, 1, 2 ... n) 

define the Bezier curve of degree n. In this way the isoscalar p art V^^^ follows the shape 
of the stiffest APR model (A18+UIX) (lAkmal et al.lll998r) up to 10 times saturation 
density uq = 0.16 fm""*. By use of a Bezier curve the model was corrected to fullfill 
the saturation point properties, like binding energy -16 MeV and compressibility Kq = 
240 MeV which were not satisfied by the orig inal A18 +UIX. With it, the most massive 
neutron star observed of about 2 Mq (Demor est et all bolO) can be created within this 
model. For the isovector part, the symmetry energy Eg, four different models sharing 
the same high density behaviour but having different slopes at saturation density have 
been introduced. The measured values of Es{no) and L (related to the slope of E.^) 
of about 30 MeV and 40-120 MeV are used for drawing the E^ curves. All of them 
satisfy the DUrca constraint that dictates that low mass neutron stars must not cool 
by DUrca process (Popov et al. 2006) related to the proton fraction x inside the star, 
which must stay below the DUrca proton fraction threshold x. Figure [T] shows these 
symmetry energy forms and the resulting proton fractions together with the DUrca 
thr eshold. Control points for y^^^ are presented in table |2] whereas for E^ can be found 
in dKubis & Alvarez-CastiU"oll2012h . 



L=40 (solid) 
L=60 (short dashed) 
L=^yO (lone dashed) 
L=I2Q (doued) 








0.12 


■\DUrca ^.^^^^ 








0.10 










0.08 










^0.06 


















0.04 
0.02 






DUrca threshold 
APR-L40-hlgh-b ■ 
APR-L60-hlgh-b ■ 
APR-L90-lilgh-b ■ 
APR-L120-hiGh-b | 


0.00 











n[fm-^] 



Figure 1 . Left. Different symmetry energy shapes for the APR-L-high-b models 
which avoid DUrca cooling for low NS masses. Right. Proton fraction as a function 
of baryon number density for these models and DUrca threshold. 



2. Neutron star properties 



To determine the neutron star properties two EoS describing its core and its crust are 
joint at equal pressure and density. For the crust the SLy EOS composed of different 
parts whose table can be found in ( IofFeli201 2) has been used. The hquid core is de- 
scribed by Eq. ([T|) and the Bezier curves described above. To de termine the crust-core 
transi tion three different methods are used which aie presented in dKubis & Alvarez-Castillo 
20121) and the resulting transition densities Uc are presented in table [T] The macro- 



scopic properties like radius, mass, moment of inertia are derived in the framework of 
General Relativity by solving the TQV equations and the crust th ickness is derived by 
the use of ndKfj) as described in ( Alvarez-Castillo & Kubisll2011 ). Figure 2 shows the 
resulting neutron star properties for each family. They all produce high enough masses 
and present thick crusts. In particular the moment of inertia carrie d by the crust do es 
not impose stringent constraints due to the glitch model restriction (iLink et al.lll999h . 
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Table 1 . Crust-core transition densities for the APR-L-high-b models. 



model 






Mc(l <r^ 2) 


APR-L40-high-b 


0.103076 


0.11012 


0.116185 


APR-L60-high-b 


0.0922071 


0.101941 


0.104926 


APR-L90-high-b 


0.0870974 


0.102523 


0.103922 


APR-L120-high-b 


0.115633 


0.142939 


0.147017 
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Figure 2. Neutron star features for the APR-L-high-b mod els. The horizont al line 
in the right lower figure represents the glitch constraint from (iLink et al.lll999h . 



It is important to mention that the inclusion of V in these models is cru- 
cial for creating massive enough neutron stars. Furthermore the particular form of 
Es fulfills the DUrca requirement. These both constraints cannot be satisfied in par- 
alell with other is oscalar parts not so stiff, hke is the case of the PAL parametriza- 
tion presented in dKubis & Alvarez-Castillol |2012|) . Therefore the models here are 
good candidates for the EoS that can be compared to microscopic approaches. Lately 
a new study that combined different laboratory measuremen ts points out values o f 
Esirio) ~ 32 MeV, L « 50 MeV with an error of a few MeV dLattimer & Limll2012h . 
From that, one may conclude that models with low L values are preferable. 
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Table 2. Bezier control points for the V isoscalar function which follows the 
APR A18+UIX EoS at high n. 



Po 




Pi 


P2 


P3 


(0.0016, 0.4649) 


(0.08, 


-20.9676) 


(0.16,-14.6553) 


(0.24, -27.2516) 


P4 




Ps 


P6 


P7 


(0.48, -29.7266) 


(0.8, 


86.6426) 


(1.12, 320.0992) 


(1.6, 1093.7393) 
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